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Radio observations of HD93129A: The earliest O star with the 

highest mass loss ? 
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Abstract. We present the results of radio continuum observations towards the open cluster Tr 14, where our main 
targets are the early-type O stars HD93129A/B and HD 93128. The observations were carried out at 3 cm (8.64 
GHz) and 6 cm (4.80 GHz) with the Australia Telescope Compact Array. Only HD93129A (type 02 If*) was 
detected; we measure flux densities of S^cm = 2.0 ± 0.2 mjy and S^cm = 4.1 ± 0.4 mjy. The resulting spectral 
index of a = —1.2 ± 0.3 (S v oc v a ) indicates predominantly non-thermal emission, suggesting HD 93129A may be 
a binary system. We propose that the observed 3 cm radio emission is mostly coming from the non-thermal wind 
collision region of a binary, and, to a lesser extent, from the thermal winds of the primary and secondary stars in 
HD93129A. At a stellar distance of 2.8 kpc, we derive a mass-loss rate M = 5.1 x 10~ 5 Mq yr _1 , assuming the 
thermal fraction of the 3 cm emission is ~0.5. 

Key words, stars: early-type — stars: individual: HD 93129A, HD 93129B, HD 93128 — stars: mass loss — stars: 
winds, outflows — radio continuum: stars 



1. Introduction 

Massive, early-type O stars at the upper end of the stellar 
main sequence are known to have strong stellar winds, 
with velocities up to ~3000 kms -1 (Prinja et al. 1990). 
Their evolution into Wolf-Rayet (WR) stars is marked by 
high mass loss rates of 10~ 7 to 10~ 5 M Q yr~ x (e.g., Lamers 
& Leitherer 1993). The transition can be described by 
three Of phases (O - 0((f)) - O(f) - Of - WR), where "f" 
stands for the presence of certain N in and He II spectral 
lines (see Mafz Apellaniz & Walborn 2002 for details; also 
Table 1). Once the H-rich outer layers of the star are lost, 
the high temperature core is exposed, and the evolution 
from Of to WR star is complete. The action of the stellar 
wind on the ambient interstellar medium creates a cavity 
of hot, rarefied gas, surrounded by a slowly expanding 
envelope (Benaglia & Cappa 1999). 

Wolf-Rayet stars have anomalously strong and broad 
emission lines, which are thought to come from the rapidly 
expanding stellar wind, but few absorption lines. With an 
average mass of 10 M©, their mass loss rate can reach 
values up to 10~ 4 Mgyr 1 . Van der Hucht et al. (2001, 
2002) cataloged well over 200 WR stars in our Galaxy, 
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more than half of which belong to the WN class, i.e. they 
show strong helium and nitrogen emission lines. A large 
fraction of the cataloged WR stars (~40%) are known to 
be binaries. 

According to Crowther et al. (1995), the most massive 
(> 60 Mq) O stars evolve from the Of stage directly into 
highly luminous late WN(+abs) and WN7 stars without 
passing through an intermediate Luminous Blue Variable 
or red supergiant phase. The slightly less massive O stars 
evolve through either of those intermediate phases into 
slightly less luminous WN8 stars (see also Langer et al. 
1994). The quantitative spectroscopic properties show a 
smooth progression from Of stars to lower excitation or 
late WN stars (WNL) (Bohannan 1990). Walborn et al. 
(1992), Lamers & Leitherer (1993) among others, sug- 
gested that transition types Ofpe/WNL, and WNL stars, 
form a sequence of increasing wind density. 

Examples of transition objects are Sk-67° 22 
(type 03IP/WN6-A, Walborn 1982), Sk-71° 34 (type 
04f/WN3, Conti & Garmany 1983), and, lately, the ones 
found in 30 Doradus (Massey & Hunter 1998). The mass 
loss rates of the transition objects are expected to lie be- 
tween those for Of and WR stars. 

Almost 400 Galactic O stars have recently been com- 
piled in the new 'Galactic O Star Catalogue' (GOS) by 
Mafz- Apellaniz & Walborn (2002). There are only about 
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a dozen stars with spectral types 03.5 or earlier; five of 
these belong to the clusters Trumpler 14 and Trumplcr 16, 
in the Carina region. Among them one star, HD93129A, 
stands out as the earliest, most luminous, and most mas- 
sive star known in our Galaxy. Walborn et al. (2002) 
re-analyzed spectra of HD93129A, previously classified 
as 03, and find that it is the prototype of a new sub- 
category and the sole representative of type 02 If* (see 
Table 1). HD93129A appears to have the most extreme 
stellar winds regarding ionization and terminal velocities 
(Walborn et al. 2002): Voc = 3200 ± 200 kms" 1 , an ef- 
fective temperature T c ff = 52, 000 ± 1000 K, and a stel- 
lar mass of 130 ± 15 M Q (Tarcsch et al. 1997). Only 
2.6" from HD93129A, Walborn (1973) discovered a sec- 
ond star, HD93129B, of type 03.5 V((f+)). Another star 
of the same type, HD 93128, is found at a projected dis- 
tance of ~ 20". The stellar parameters for all three stars 
arc listed in Table 1. 

In the radio regime, free-free thermal emission from 
the ionized wind of early-type OB and Wolf-Rayet stars 
can be modeled (Wright & Barlow 1975, Panagia & Fclli 
1975) yielding a spectral index (a) of typically ~ 0.6 — 0.7. 
Generally, the radio flux densities are of the order of mJy 
for stars at distances up to 2 kpc. About 80 WR stars 
(Abbott et al. 1986, Chapman et al. 1999) and -120 O 
stars (Bieging et al. 1989, Scuderi et al. 1998, Benaglia et 
al. 2001b, etc.) have been observed. Roughly half of them 
were detected at least at one frequency. When spectral 
indices could be measured, ~40% showed non-thermal or 
composite indices. 

In recent years the hypothesis was proposed (Chapman 
et al. 1999, Dougherty & Williams 2000) that non-thermal 
emission could be indicative of binarity. In that case, 
part of the emission comes from the colliding wind re- 
gion (CWR) of two early massive stars. Chapman et al. 
(1999) observed 36 WR stars with the ATCA at 3, 6, 13 
and 20 cm, obtaining an angular resolution of ~ 1" at 
3 cm. Dougherty & Williams (2000) discussed the spectral 
index behaviour of 25 WR stars, by taking into account 
observations carried out mainly with the VLA, ATCA, 
MERLIN and the WSRT. Benaglia & Romero (2003) pre- 
sented a table with the characteristics of the non-thermal 
or composite WR emitters. 

Unfortunately, in some particular cases where the non- 
thermal emission is strong, not enough spectroscopic data 
are currently available to investigate the binary status. For 
example, CD-47° 4551 (CPD-47 2963), an 04 111(f) star 
(Walborn 1982), was observed with the ATCA (Benaglia 
et al. 2001b) and detected with flux densities of 1.77±0.05 
mJy at 3 cm and 2.98 ± 0.05 mJy at 6 cm (a = -0.89 ± 
0.06). If the star is a binary with another WR or OB star, 
the non-thermal emission could be produced at the CWR 
of both stars. 

Non-thermal radio emission can also be produced by 
the stellar wind of a single star, as described by White 
(1985). In his model, electrons are accelerated to relativis- 
tic energies by shocks in the wind near the star, and emit 
radio radiation. For strong shocks, the expected spectral 



index is +0.5. Though this is similar to the thermal index 
value, the non-thermal flux density coming from a single 
stellar wind will suffer a break at a frequency of a few 
GHz, and behave as S(v) oc log(^) for larger frequencies, 
while the thermal flux will obey a a = 0.6 — 0.7 law up to 
the infrared range. 

Radio observations are very useful to derive the mass 
loss rates of early- type stars (Leitherer et al. 1995). The 
source HD93129A represents one of the most interesting 
candidates because it is the only known example of the 
hottest supergiant class. With the help of high resolution 
radio data, properties of the wind can be unveiled. In de- 
tecting the stellar wind, we intend to identify its radio 
spectrum, look for a possible colliding wind region, and 
compare the mass loss rate with those derived from other 
methods. The observations will also allow us the study of 
other early-type stars in the region. 

The paper is organized as follows: Section 2 summa- 
rizes the characteristics of the target sources relevant to 
this study. Section 3 describes the observations and data 
reduction. In Section 4 we present the results, followed by 
the discussion in Section 5. Section 6 contains our conclu- 
sions. 

2. The target stars in Trumpler 14 

The field of view of the present data (see Section 3) corre- 
sponds to the cluster Trumplcr 14. Lists of the catalogued 
members can be found in Massey & Johnson (1993), and 
recently DeGioia-Eastwood et al. (2001). A wide range 
of spectroscopic or photometric distances can be found 
for Trl4 in the literature. Humphreys & McElroy (1984) 
derived a distance of 3.47 kpc. Morrell et al. (1988) dis- 
cussed two possible distances, 2.8 (variable R, from Tapia 
et al. 1988) and 3.45 kpc (R = 3.2), depending on the 
value taken for the ratio of total to selective absorption 
R. Massey & Johnson (1993) and Vazquez et al. (1996) 
found similar values of 3.2 and 3.1 kpc, respectively. Most 
recently, Tapia et al. (2003) computed spectroscopic paral- 
laxes to derive a distance of 2.8 kpc. We adopt the latter 
value throughout this paper. At this distance, the pro- 
jected cluster size of 140" (Tapia et al. 2003) corresponds 
to 1.9 pc. Chandra observations towards the rj Carina re- 
gion (Evans et al. 2003) revealed strong X-ray emission 
from the stars HD 93129A/B. Their combined X-ray lumi- 
nosity is L x (0.5 - 2.4keV) = 2.9 x 10 33 ergs" 1 , adjusted 
to a distance of 2.8 kpc. 

Penny et al. (1993) presented radial velocity mea- 
surements for seven stars in Tr 14, with mean values of 
-4.4 ± 7.4 kins" 1 for HD93129A, 7.2 ± 2.4 kms" 1 for 
HD 93129B and 7.3 ± 3.3 km s" 1 for HD 93128. They con- 
cluded that none of them is a strong binary candidate 
(see also Conti et al. 1977, 1979). Based on the standard 
deviation in the mean stellar radial velocities, they esti- 
mated that the typical random motions in this cluster are 
8.5 kms -1 . For HD93129A Penny et al. placed an upper 
limit of ^10 kms -1 on the semiamplitude of any short- 
period binary. To search for possible long-period bina- 
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Table 1. Stellar parameters adopted for the target stars 





HD 93129A 


HD 93129B 


HD 93128 




a, 8 (J2000) 


10 43 57.462 (a) 


10 43 57.638 (b) 


10 43 54.372 (c) 


h m s 




-59 32 51.27 


-59 32 53.50 


-59 32 57.37 


/ 51 


Spectral Class. 


02 If* (d) 


03.5 V((f+)) (d) 


03.5 V((f+)) (d) 




Voc 


3200 ± 200 (o) 


3070 (f) 


3070 (f) 


kms" 1 




52000 ± 1000 (e) 


52000 ± 1500 (g) 


52000 ± 1500 (h) 


K 


log (L/L ) 


6.4±0.1 (o) 


5.6 (g) 


5.7 (s) 





References: (a): Hog et al. 2000; (b): Worley & Douglass 1997; (c): see Mafz Apellaniz & Walborn 2002; (d): Walborn et al. 
2002 (f* means the Niv 4058A emission lines are stronger than the Niii 4640A emission lines, and Hell 4686A is observed in 
emission, while ((f+)) means weak N ill 4634-4640-4642A emission lines, strong He n 4686A absorption lines, plus Si IV 4089- 
4116A emission lines); (e): Taresch et al. 1997; (f): interpolated from Prinja et al. 1990; (g): see text; (h): Vacca et al. 1996. 



ries they obtained speckle observations, with negative re- 
sults. Surprisingly, recent HST observations indicate that 
HD93129A is a binary (Walborn 2002) with a separation 
of 55 mas (Walborn 2003, priv. comm.; Nelan et al., in 
prep.), which could have been detected by Penny et al. 
(1993). 

The optical positions of the three earliest stars in the 
observed field were taken from the GOS catalog, and 
come from three different sources, varying in precision: 
the Tycho-2 catalog (Hog et al. 2000) lists a very accu- 
rate position (a <~ 5 mas) for HD93129A, whereas the 
position of HD93129B originates from the Washington 
Double Star catalog (Worley & Douglass 1997) and that 
of HD 93128 was measured from NTT data, according to 
Mafz- Apellaniz & Walborn (2002). For a summary of the 
stellar parameters see Table 1. 

HD 93129A. The spectrum of this early 02 If* star is 
very similar to Wolf-Rayet stars of spectral type WN6 or 
WN7 (Walborn et al. 2002). Taresch et al. (1997) stud- 
ied the spectrum from optical to FUV range, concluding 
that its strong wind is radiatively driven. By fitting an 
Ha profile they obtained a mass loss rate of 1.8 x 10~ 5 
M©yr -1 compared to 2.2 x 10~ 5 M Q yr _1 derived from 
Ha measurements by Puis et al. (1996; assuming T e g = 
50,500 K). 

HD93129B and HD 93128. The effective tempera- 
ture of HD 93128 is 52000 ± 1500 K (Vacca et al. 1996). 
For HD 93129B, which has the same spectral type as 
HD 93128, we adopt the same value. The stellar luminosi- 
ties were derived taking the absolute visual magnitudes 
from Walborn et al. (2002). 

HD 93250. This star is a member of the Tr 16 cluster, 
located ~ 10' East of HD 93129A/B (well outside our pri- 
mary beam). Leitherer et al. (1995) reported a 3cm flux 
density of 1.36 ± 0.17 mJy and a 3er upper limit of 3.57 
mjy at 6 cm. Assuming a purely thermal spectrum, a dis- 
tance of 2.2 kpc, a wind velocity of = 3200 kms -1 and 
an effective temperature of T cff = 50,000 K (T e = 0.4 T cS 



= 20,000 K according to Drew 1989) Leitherer et al. de- 
rived a mass loss rate of 4.1 x 10~ 5 M Q yr _1 . Considering 
a re-classified spectral type of 03.5 V((f+)) (Walborn et 
al. 2002) and a distance of 2.5 kpc to Trl6 (Tapia et al. 
2003) we revise the mass loss rate to 5 x 10~ 5 MQyr" 1 
(see Section 4.1). 



3. Observations and Data Reduction 

The observations presented here were obtained with the 
Australia Telescope Compact Array (ATCA) in January 
2003 using the 6B configuration. The target field, which 
contains three massive Of stars, was observed simultane- 
ously at two frequencies: 8.64 GHz (3 cm) and 4.80 GHz 
(6 cm), with a total bandwidth each of 128 MHz over 32 
channels. The maximum baseline was 6 km, resulting in 
angular resolutions around 1" —2" . The size of the primary 
beam is 10' at 6 cm and 5' at 3 cm. The pointing position 
was a,<5(J2000) = 10 h 43 m 56 s , -59° 32' 51", slightly off- 
set from the three stars HD 93129A/B and HD 93128. The 
field was observed during intervals of 8 to 15 min, depend- 
ing on the atmospheric phase stability, interleaved with 
short observations of the nearby phase calibrator 1045- 
62, and spanning an LST range of 12 h. The total time on 
the target field was nearly 300 minutes. 

The data were reduced in miriad using standard pro- 
cedures. The flux density scale was calibrated with obser- 
vations of the primary calibrator PKS B1934-638, assum- 
ing flux densities of 2.84 Jy at 3 cm and 5.83 Jy at 6 cm 
(Reynolds 1994). After data calibration the visibilities 
were Fourier-transformed using 'natural'- and 'uniform'- 
weighting. To overcome confusion due to diffuse emission 
from extended sources in or towards the observed fields, 
the shortest baselines were removed. For the same rea- 
son 'uniform'-weighting resulted in the lowest noise levels: 
0.11 and 0.21 mJybcam -1 at 3 and 6 cm, respectively. The 
synthesized beam is 1.1" x 1.0" at 3 cm and 1.9" x 1.7" at 
6 cm. 



Paula Bonaglia & Barbel Koribalski: The HD 93129A field 









1 1 1 1 1 1 £ ~ u ' 1 1 1 1 




30" 






° . ' p ' ' 4 


o n , ■ 


40" 




1 o 





<- ! 


50" 


- 


i 


^- HD93129A ° 

TjnAn 1 ODD ^ G 






? 


33'00" 






o 




HD93128 . tr 




59 - 33'l0" 


- • 









o 



10°43 m 59 s 58 s 57 s 56 s 55 s 54 s 53 s 

RA (J2000) 

Fig. 1. ATCA 3 cm radio continuum image towards Tr 14. 
The optical positions of three stars (see Table 1) are 
marked. The contour levels arc -0.22, 0.22 (2a), 0.44, 
0.66, 0.88, and 1.10 mJybeam -1 . The synthesized beam 
(1.1" x 1.0") is displayed at the bottom left. 
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Fig. 2. ATCA 6 cm radio continuum image towards Tr 14. 
The optical positions of three stars (see Table 1) are 
marked. The contour levels are -0.42, 0.42 (2a), 0.84, 
1.26, 1.68, and 2.10 mJybeam -1 . The synthesized beam 
(1.9" x 1.7") is displayed at the bottom left. 



Table 2. Flux densities (or 3a upper limits) of the target 
stars. 



Star 


5*301X1 


Sqctci 




(mjy) 


(mjy) 


HD 93129A 


2.0 ±0.2 


4.1 ±0.4 


HD 93129B 


< 0.33 


< 0.63 


HD 93128 


< 0.33 


< 0.63 



4. Results 

Figs. 1 and 2 show the central part of the resulting 3 and 
6 cm radio continuum images. The optical positions of the 
three target stars (see Table 1) are marked. 

Only one radio source, HD 93129A, is detected. Using 
the miriad task imfit we fit the position and flux den- 
sity (of a point source) at both frequencies resulting 
in 'a,(5(J2000) = 10 h 43 m 57 s .45, -59° 32' 51.36" , 2.0 ± 
0.2 mjy at 3 cm, and a,<5(J2000) = 10 h 43 m 57 s .47, 
-59° 32' 51.23" , 4.1±0.4 mjy at 6 cm. From the 3 and 6 cm 
flux densities we derive a spectral index of a = —1.2 ± 0.3 
(S v oc v a ) indicating predominantly non-thermal emis- 
sion. The measured flux densities for HD 93 129 A and the 
3cr upper limits for HD 93129B and HD 93128 are listed in 
Table 2. 

The separation between the fitted 3 and 6 cm radio 
continuum positions for HD 93 129 A is 0.2", slightly larger 
than the positional uncertainty (which is approximately 
the size of the synthesized beam divided by signal-to-noise 
ratio of the source). The radio and optical positions agree 
within the given uncertainties. 

The Simbad database lists 250 Galactic objects within 
a radius of 2.5' around our pointing position (see 



Section 3), at least 11 of which are O-type stars, includ- 
ing the three target stars. Table 3 lists the eight additional 
stars, none of which was detected. The same 3a upper lim- 
its apply. 

4.1. Mass loss rates 

According to Wright & Barlow (1975), the stellar mass 
loss rate, M, of an optically thick ionized wind can be 
expressed in terms of the measured thermal flux density, 
S v , at a frequency v, from radio to IR ranges, as: 

C 3 / 4 J3/2 

M = 5.32 x 10- 4 bv Moyr- 1 . (1) 

Here d is the stellar distance in kpc, the wind termi- 
nal velocity in kms -1 , S v the thermal flux density in mjy, 
and v the frequency in Hz; \i stands for the mean molec- 
ular weight of the ions, 7 is the mean number of electrons 
per ion and Z the r.m.s. ionic charge. For a fully ionized 
H-Hc gas, [i — 1.38. In a first approximation, Z w 7 ~ 1. 
We will take those values for all calculations throughout 
the text. The Gaunt factor <?„ can be approximated by: 

* = 9.77(l + 0.13 ,„ 6 M^kff) . p, 

From the derived spectral index of a = — 1.2 ± 0.3 
we know that the observed radio continuum emission 
consists of thermal and non-thermal components. Their 
relative contributions (and locations) are unknown 
since the detected source is essentially unresolved. If 
/t is the fraction of the thermal emission at 3 cm with 
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Name 




a, S (J2000) 


Spectral type 


M cxp 


"-'expjScm 






(h m s) 


(°'") 




M yr" 


-l 


mjy 


CI* Trumpler 14 MJ 165 


10 43 55.43 


-59 32 49.8 


08 V (a) 


8.5 x 10" 


-7 


0.010 


CPD-58 2620 




10 43 59.92 


-59 32 25.4 


06.5 V((f)) (b) 


1.7 x 10" 


■6 


0.016 


CI* Trumpler 


14 MJ 150 


10 43 53.75 


-59 33 29.0 


09: V (a) 


4.5 x 10" 


-7 


0.005 


CI* Trumpler 


14 MJ 127 


10 43 48.82 


-59 33 24.8 


09 V (a) 


4.5 x 10" 


-7 


0.005 


CPD-58 2611 




10 43 46.69 


-59 32 54.8 


06 V((f)) (b) 


2.1 x 10" 


■6 


0.021 


HD 93160 




10 44 07.27 


-59 34 30.8 


06 III((f))( b ) 


4.9 x 10" 


■6 


0.064 


HD 93161AB 




10 44 08.90 


-59 34 34.9 


06.5 V((f)) (b) 


1.7 x 10" 


■6 


0.016 


LS 1832 




10 44 06 


-59 35 09.0 


06 (c) 


2.1 x 10" 


■6 


0.021 



References: (a): Massey & Johnson 1993; (b): Mai'z-Apellaniz & Walborn 2002; (c): Reed & Beatty 1995. 



Sjcm = /t x ^cm, /t = — 1, then the mass loss rate 
of HD93129A is M — (/ T ) 3/4 x 8.6 x 10~ 5 Moyr" 1 . 
For /t ^ 1 this would be one of the highest mass loss 
rates ever derived for a massive, early-type star. For 
HD93129B and HD 93128 we derive 3cr upper limits of 
M < (/t) 3/4 x 2.13 x 1(T 5 M yr- 1 . 

Given a stellar mass, luminosity, effective tempera- 
ture, and wind terminal velocity, the models of Vink et 
al. (2000) can be used to compute a theoretical — or ex- 
pected — mass loss rate, M oxp . Using Vink's program 
(see astro.imperial.ac.uk/~jvink/) and the values given 
in Table 1, we calculate M cxp = 2.12 x 10~ 5 Meyr" 1 for 
HD93129A. While this value is in good agreement with 
the Ha-derived mass loss rates (Puis et al. 1996, Taresch et 
al. 1997), it is a factor four smaller than that derived from 
the measured 3 cm radio continuum flux density, suggest- 
ing / T = 0.15. For HD 93129B and HD 93128 we calculate 
M exp of 0.3 x 10~ 5 M/gyr- 1 . 

Table 3 lists the expected stellar mass loss rates for 
the other eight O stars in the field, estimated as for the 
earliest stars and assuming a distance of 2.8 kpc. The lu- 
minosities, masses, and effective temperatures were taken 
from Vacca et al. (1996). The terminal velocities were in- 
terpolated from Table 3 of Prinja et al. (1990) according to 
spectral type. The expected thermal flux density at 3 cm 
was computed inverting Eq. (1), and is listed in the last 
column. 



5. Discussion 

5.1. A CWR between HD 93129 A and B ? 

The distance of D — 2.6" = 7300 AU between coordinates 
of HD 93 129 A and B is a lower limit to their absolute 
separation. If a colliding wind region exists between the 
two stars, it would be located at projected distances r A , 
r B (see Eichler & Usov, 1993): 



The parameter i] is defined in terms of the wind ter- 
minal velocities and the stellar mass loss rates M: 



'/ — 1 ■ 

M A X Voo^a 

Using the derived theoretical values for M A and Mb, 
we find r\ = 0.14. Thus r A = 1.9" w 5350 AU, and r B 
0.7" w 1950 AU. The putative CWR would be centered 
around a, <5(J2000) C wr = 10 h 43 m 57.59 s , -59° 32' 52.9" , 
but no emission was detected with the ATCA. The same 
3<t upper limits apply (see Table 2). 

Assuming a spectral index of a = —0.5 (Eichler & 
Usov 1993), we derive an upper limit to the non-thermal 
radio luminosity from a CWR between stars A and B of 
L b s < l-4x 10 28 ergs -1 . The expected non-thermal lumi- 
nosity from a possible A/B binary is difficult to estimate 
due to the poor knowledge of system parameters such as 
size of the CWR and the magnetic field. Moreover, ac- 
cording to Wright & Barlow (1975), the radii of the radio 
photospheres at 3 and 6 cm would be ~30 and 60 AU, 
(<< r A , r B ), respectively, if M = 5 x 10~ 5 M yr" 1 . The 
existence of a CWR between HD 93129A and B is very 
unlikely (see Section 5.4). 

5.2. A CWR between the binary components of 
HD 93129 A? 

Walborn (2002) recently reported that HD93129A is a 
binary system. Using the HST Fine Guidance Sensor 
Interferometer, Nelan et al. (2003, in prep.) resolved it into 
a 55 mas binary with a magnitude difference of 0.9 mag 
(Walborn 2003, priv. comm.). In the binary system, we 
will call "a" the 02 primary, and "b" , the secondary. The 
total V magnitude measured for HD93129A is 7.2 mag 
(Hog et al. 1998). Given that U a = 7.6 mag, and V h = 8.5 
mag (Walborn 2003, priv. comm.), the secondary has a 
magnitude similar to that of HD93129B (V = 8.9 mag) 
and HD 93128 (V = 8.8 mag) (Walborn et al. 2002). The 
four stars are probable coeval, and their spectral types 
should be similar (E. Nelan et al., in preparation). 

The binary separation (55 mas = 154 AU) implies that 
the stellar winds are likely to interact. In fact, there is 
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a high probability that the system is an O+O colliding 
wind binary. Additional support to this idea can be found 
by studying the radio spectral index (see Section 5.3). 
Assuming that the companion is an 03.5 V star with 
the same terminal velocity, effective temperature, lumi- 
nosity, and expected mass loss rate than HD93129B (see 
Table 1), we find n = 0.14, r Aa = 112 AU, and r Ah = 42 
AU. The binary separation of HD 93129A is similar to that 
of WR146 (210 AU, Setia Gunawan et al. 2000), thus it 
is reasonable to assume a similar extent for the CWR, i.e. 
~50 AU. We estimate the equipartition magnetic field for 
the CWR according to Miley (1980) and find B C wr ~ 4 
mG. Using the magnetic field structure given by Eichler & 
Usov (1993), the surface magnetic field of HD 93129Aa can 
be extrapolated to -B* ~ 100 G, if the stellar rotational 
velocity is O.luoo (Conti & Ebbets 1997). This estimate 
is consistent with the few observational values available 
(Donati et al. 2001, 2002). 

Substituting the parameters of HD 93129Aa and Ab 
into Eqs. (15) and (16) from Eichler & Usov (1993) we 
obtain: 

L rad w 7.4 x 10 31 (j^) (i^r) ergs- 1 (5) 

for the expected non-thermal radio luminosity of the 
CWR, and £ w 5 x 10~ 5 . Taking typical values of /3 ~ 0.1, 
e <~ 10~ 2 , we calculate L rad « 4 x 10 29 ergs -1 . 

From the observed flux densities we obtain a radio lu- 
minosity of L obs = 1.12 x 10 29 ergs" 1 for HD93129A be- 
tween 4.8 and 8.64 GHz, which is in good agreement with 
the total expected one. 

In a colliding wind binary system like the one described 
here, a total X-ray luminosity due to Bremsstrahlung can 
be computed. Following Usov (1992) we obtain 2.2 x 10 33 
ergs -1 , very close to the observed X-ray luminosity by 
Evans et al. (2003). 

5.3. Spectral index considerations 

The spectral index of —1.2 ±0.3 derived from the observa- 
tions presented here is much steeper than that estimated 
by White (1985), a — +0.5, for synchrotron emission from 
a single star. Steeper indices can result in the CW scenario 
as the effect of the softening of the particle spectrum due 
to inverse Compton (IC) losses (Benaglia & Romero 2003). 
At a CWR the accelerated electrons follow a power-law 
spectrum N(E) oc E p , with p = —2. Both synchrotron 
and IC processes help electrons to lose energy. At high 
energies, the radiative cooling would be dominated by IC 
losses, due to the strong photon fields at the CWR. At the 
energy E h at which the cooling and escape times are equal 
(e.g. Longair 1997, p. 281), IC losses produce a break in 
the spectrum, from an index p to p' = p— 1 . The same pop- 
ulation of rclativistic electrons is capable of losing energy 
via synchrotron emission, producing a flux S(v) oc v a , 
where a = (p+l)/2. Thus, for E < E hl a w —0.5, and for 
E > £b, at w — 1. This last value seems to be consistent 
with the spectral index derived from the present observa- 



tions. The break energy can be computed if the follow- 
ing set of parameters is known (see Benaglia & Romero 
2003): the mass loss rates, wind terminal velocities, and 
luminosities of both components, and the magnetic field 
at the CWR, together with the size s of the non-thermal 
source. Replacing the values adopted or obtained above, 
the break frequency in the radio range will be a few GHz. 
ATCA observations at 1.4 and 2.4 GHz are necessary to 
test this scenario. 

The relative importance of the processes already men- 
tioned can be measured by estimating the time needed 
for an electron to lose all its energy, either in the form 
of synchrotron radiation (t syn ) or IC scattering (tic)- The 
importance of synchrotron and inverse Compton radiation 
for electrons of the same energy can be measured comput- 
ing the ratio between the corresponding lifetimes (Longair 
1997, p.274): 

tsyn ^rad 

he U mag 

where £/ ra d and C/ mag are the energy densities of the radia- 
tion field and of the magnetic field respectively. Using the 
observed radio luminosity L bs, we get U la d — 2.5 x 10 -9 
cV uT 3 . The variable J7 mag is 6.4x 10~ 15 cV m~ 3 . The cor- 
responding time ratio is about ~ 4 x 10 5 , thus t syn » tic- 
Following this last estimate, and under the conditions 
mentioned above, IC losses should dominate over syn- 
chrotron ones. 

5.4. Colliding wind emission from other binaries 

There are only a few examples (WR 146, WR 147, Cyg 
OB2 No. 5) of early-type binary systems in the litera- 
ture where radio emission from the CWR has been ob- 
served. WR 146 is a WC6+08 binary (van der Hucht 
2001). Dougherty et al. (2000) observed the stars with 
MERLIN at 5 GHz twice, and with the VLA at 1.4, 5, 8.5 
and 22 GHz. The combination of their MERLIN data with 
the 22 GHz VLA data allowed them to detect two ther- 
mal sources, coincident with the WR and OB stars, plus 
a non-thermal one in-between, nearer to the OB star, and 
possibly coming from the CWR. At a distance of 1.25 kpc 
(van der Hucht et al. 2001), the separation of the stellar 
components is ~210 AU, and the size of the non-thermal 
source is ~50 AU (Setia Gunawan et al. 2000). The flux 
densities quoted in the literature are: 31.4+ 0.4 mJy and 
28.5 ± 0.3 mJy at 5 GHz for the source identified as the 
CWR, 7.0 ±1.3 mJy at 22 GHz for the emission possi- 
bly related to the WR star, and 10.4 ± 1.0 mJy from the 
OB star. The ratios of the non-thermal to thermal flux 
densities are ~3 for the O star, and ^4 for the WR star. 

Contreras & Rodriguez (1999) presented 8.5 GHz- VLA 
observations of the star WR 147, a binary system formed 
by WN8 (h)+B0.5 V components (van der Hucht et al. 
2001), at a distance of 650 pc (van der Hucht 2001). They 
detected a southern, thermal radio source coincident with 
the WR component, and a northern, non-thermal source. 
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The latter is interpreted as originating at the inter; 
zone between the winds of the WR and OB stars, 
averaged measured flux densities are 28.4 ± 0.5 mJy 
the thermal source, and 10.4 ± 0.5 mJy from the 
thermal source. The stellar separation is ~420 AU 
Gunawan et al. 2001). The non-thermal to therma 
density ratio is ^0.4. 

Cyg OB2 No. 5 seems to be formed by three 
an 07 Ia+Ofpe/WN9 contact binary, and a B0 V s 
a larger orbit (Contreras et al. 1997, Rauw et al. '. 
located at 1.8 kpc (e.g. Waldron et al. 1998). The sc 
tion between the contact pair and the third compon 
about 1700 AU. Radio observations by Persi et al. ( 
and Miralles et al. (1994) at 1.4, 5, 8.5 and 15 GHz, si 
that the source coincident with the contact pair pres 
flux variations, defining a "high" emission state (S%. 5< 
5-7 mJy, and a ~ 0), and a "low" emission state (Si 
= 1-2 mJy, and a <~ 0.6). Contreras et al. (1996, 
performed VLA observations at 5 and 8.5 GHz, det< 
a non-thermal source (a = —2.4 ±0.6) near the B0 \ 
and proposed it was a CWR in the system. They mea 
a flux density at 8.5 GHz of 0.31 ± 0.02 mJy. The r£ 
the emission from the weaker non-thermal source t 
emission of the contact pair at the lower state, is 
0.2. 

The system HD93129A presents a binary sepa] 
similar to that of WR146. The ratio of non-thern 
thermal emission is unknown. If the magnetic field ; 
CWR is equal to the equipartition field derived abo\ 
mG, then it is similar to those derived for WR 147 (£ 
and WR 146 (25 mG); see Benaglia & Romero (2003; 
scarcity of the examples and the huge differences i 
ratios of CWR-non-thermal to stellar-thermal flux i 
ties measured for the three cases from the literature, 
the results of the comparison inconclusive. 

5.5. The radio spectrum of HD93129A 

To separate the thermal and non-thermal contribi 
to the spectrum of HD 93129A, further radio observ: 
are needed at lower and, particularly, at higher radio rre- 
quencies. These can be obtained with the ATCA at 1.4, 
2.3, 16-25, and 85-91 GHz. The system can be resolved 
with the Long Baseline Array, but the flux densities of 
the thermal winds are likely too low for a detection, and 
an absolute position for the CWR would be difficult to 
obtain. 
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Fig. 3. a) Measured and predicted radio continuum flux 
densities for HD 93 129 A in the frequency range from ~1 to 
25 GHz. The two measurements at 4.8 and 8.64 GHz are 
indicated by black squares. Assuming that the flux den- 
sity measured at 8.64 GHz, 5*301117 consists to equal parts 
of thermal and non-thermal emission (/t = 0.5), and that 
the thermal emission has a spectral index of cut = 0.6, we 
can fit the observed flux densities by adding a non-thermal 
component with cunt = — 2.1±0.3. The dotted and dashed 
lines show the thermal and non-thermal radio emission, 
respectively. The solid line represents the combined ther- 
mal and non-thermal radio emission for this model, b) As 
Fig. 3a, but assuming /t = 0.15. In this case, we can fit 
the observed flux densities by adding a non-thermal com- 
ponent with afjT = — 1.4 ± 0.3. 



If the actual mass loss rate of HD93129A is M — 
2.1 x 10~ 5 M Q yr" 1 , i.e. f T = 0.15 (sec Section 4.1), 

the thermal and non-thermal contributions to the 3 cm In contrast, for /t = 0.5 we find «nt ~ —2.1 ± 0.3, 

radio flux density are ~0.3 and <~1.7 mJy, respectively, similar to that derived for other Of systems (e.g., Cyg OB2 

Assuming a spectral index of qt = 0.6 for the thermal No. 5, Contreras et al. 1996, 1997). At 16, 25 and 90 GHz 

component, we find the spectral index for the non-thermal we expect <~1.7, <~2, and ^4.1 mJy (see Fig. 3b). In this 

component is «nt ~ —1-4 ± 0.3. The expected total flux case, the radio derived mass loss rate of the system is 5.1 x 

densities (see Fig. 3a) arc <~1.2, <~ 1, and <~1.3 mJy at 16, 10~ 5 M Q yr" 1 . The two scenarios can be distinguished 

25 and 90 GHz, respectively. with high-frequency ATCA synthesis observations. 
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6. Conclusions 

The field of view towards the Galactic cluster Tr 14 has 
been observed with the ATCA in the radio continuum at 
4.80 GHz (3cm) and 8.64 GHz (6cm). Only one source, 
the 02 If* star HD 93129A, was detected with flux densi- 
ties of 2.0 ± 0.2 mJy at 3 cm, and 4.1 ± 0.4 mJy at 6 cm, 
resulting in a very steep spectral index of a — —1.2 ± 0.3. 
This star, recently reported as binary, is the earliest known 
O-star and potentially the most massive system in the 
Galaxy, followed by rj Carina. The observed radio emis- 
sion is most likely coming from (a) the thermal winds of 
both binary components and (b) the non-thermal emission 
from the colliding wind region. The negative spectral in- 
dex of HD 93 129 A suggests that the latter is the dominant 
emission mechanism. 

We derive a mass loss rate of (/t) 3 ^ 4 x 8.6 x 10~ 5 M Q 
yr^ 1 for HD 93129A and 3cr upper limits of (/ T ) 3/4 x 2.13 x 
10~ 5 M Q yr- 1 for HD93129B and HD 93128, where / T is 
the thermal fraction of the 3 cm radio continuum emission. 

Although the secondary star in the binary HD 93 129 A 
remains yet unclassified, the presence of intense non- 
thermal emission and high radio luminosity together sug- 
gest it is another massive, early-type star. 

The lack of radio emission near HD93129B, and the 
extremely large separation of HD93129A/B (7300 AU at 
2.8 kpc) conspires against the existence of a colliding wind 
region in-between the stars. 

It is crucial to perform observations (i) at lower fre- 
quencies, where the thermal contributions should be neg- 
ligible, and an accurate non-thermal index can be derived; 
this would help to investigate the processes responsible 
for the radiation detected; and (ii) at high frequencies, at 
which non-thermal emission decreases, and the estimation 
of the thermal flux will enable us to derive the value of 
the true mass loss rate of HD 93129A. 
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